GAYLUSSACIA 






PHYLOGENETIC ANALYSES OF THE 


COMPLEX (ERICACEAE: VACCINIEAE) 
AND MORPHOLOGICAL CHARACTERS 




Michael T. Gajdeczka, Kurt M. Neubig 

Department of Biology 
220 Bartram Hall 
P.O. Box 118526 
University of Florida 
Gainesville, Florida 32611, USA. 


michcio@ufl.edu; kneubig@flmnh.ufl.edu; w. 



er S. Judd 



W. Mark Whitten, Norris H. Williams, Kent D. Perkins 

University of Florida Herbarium 
379 Dickinson Hall 


P.O. Box 110575 
Florida Museum of Natural History 

University of Florida 
Gainesville, Florida 32611, U.5.A. 
ufl.edu; orchid@flmnh.ufl.edu; 




ABSTRACT 

The three closely related species of blue huckleberries, i.e. Guv/ussacuz fwndosa , G. tomentosa, and G. nana, which make up the Gaylus- 
wauJrondosa complex (sect. Decamerium , subsect. Frondosae'). have been variously circumscribed and classified by North American 
siemausts. This study aims to resolve the species relationships of the Gaylussacia frondosa complex through phylogenetic analyses of 
morphological characters and nuclear ribosomal internal transcribed spacer (nrlTS) sequences. Gavlussutiu nana forms a distinct clade 
with support in all analyses. Contrary to some earlier treatments, our results strongly s upport distinguishing G. nana from G. towientosu, 
though a few putative hybrids between these species were discovered. The monoph ly of G. tomentosa is not strongly supported, but 
ts upheld in the morphology-based analysis and the combined (nrlTS + morphology) analysis, and a G.jfftHuioStt + G. tnrnentosu clade is 
strongly supported in the nrlTS and combined analyses. In contrast, the morphological topology groups {?«frondosa with G. narUL None 
analyses support a G. Jrondosa clade, but the specimens representing G. /roudosu do form a clade in some of the most parsimom- 
i-^5 m the combined analysis, while in others they constitute a paraphyletic complex (in a clade with G. tomentosa). Chromosome 
unts y previous workers indicate that G. tomentosa is teiraploid, whereas G. frondo$& and G. nana are diploid. Gtivlussuciti/rondosa (in 
Erast to C. tomentosa and G. nana) tends to be a taller plant and has a more northern distribution. Based on this evidence, we sugge 
51 three taxa be considered distinct species (each likely representing a separate evolutionary lineage). 

DS * ssaciaj rondosa, G. tomentosa , G nana, Ericaceae, Gaylussacia sect. Decam 
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r ' es P ec tcs intimamente relacionadas de arandaneras, i.e., Gaylussacia frondosa, G. tomenfiisa, y G. nana, que constituyen 

fi&Qdosit (sect. Decamerium, subsect. Frondosae), han sido circunscritas y clasificadasvariadamente por los sistcmaticos norte 


nos. Este estudio trata de resolver las relaciones entre las especies del complejo Gaylussacia fwndosa mcdiante analisis filogeneticos 
ffes morfologicos y secuencias del espaciador nuclear ribosomico interne (nrlTS). Gaylussacia nana forma un clado distmto con 


^ en todos los analisis. Contrariamente a algunos tratamientos previos, nuestros resullados soportan fuertemente la distinct>n 

l ^^ towentosa, aunque se encontraron algunos hibridos putativos entre estas especies. La monofilia de G tomentosa no esta 

m ente sonnrt^^^ ~- -- « « — • * 1 J / — , * rr morfologia), y un clado G. frondosa + G. 



sino que esta apoyada en el analisis morfolbgico y el combmado (nrlTS + 

_ _ es * a ^ uer temente soportado en los analisis de nrlTS y los combinados En contraste, la topologia morfolbgica agrupa G. frondosa 
U , hJinguno de los anahsis soporta un clado de G. frondosa, pero los especimenes que representan G. fwndosa forma un clado 

^ 0s wholes mas parsimoniosos del analisis combinado, mientras que en otros conslituye en complejo parafiletico ten un 

iCc *conG f/ini i. | 

Bl _, nuosa). Los recuentos cromosomaticos realiz 

"films que G, frondosa y C. 

^ytieneuna distribuc 

na re prcsentando probablemente una linea evolutiva separada) 



s por investigaaores anteriores mdican que G. tomentosa es tetraploide, 
nana son diploides. Gaylussacia frondosa (en contraste con G. tomentosa y G. nana) tiende a ser una planta 
ion mas septentrional. Basados en estas pruebas, sugerimos que los tres taxa sean considerados especies 
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Gaylussacia Kunth is a New World genus of 53 species of mainly understory shrubs occurring in raesic to 
xeric woodlands and shrublands, as well as acidic bogs. The genus is differentiated trom the closely-related 
Vaccinium L. by ten-locular ovaries, drupaceous fruits containing ten pits, the presence bf resin glands on the 
leaves (except for G. brachyccra (Michx.) A. Gray), and the lack of staminal spurs (Duncan & Brittain 1966; 
Luteyn et al. 1996; Palser 1961; Wood 1961). However, the phylogenetic relationships between membersof 
the two genera are still uncertain i Kron et al. 2002b). Vocciflium was shown to be non-monophyletic as tradi¬ 
tionally circumscribed, with Gaylussacia nested within it (Kron et al. 2002a). These results are corroborated 
by the ambiguous generic identity of G. brachyccra, which, in addition to being eglandular and similar in 
form to Vdcci mum species, has been shown to be sister to the rest of Gaylussacia (Floyd 2002). Finally, some 
Vacciriium species have ovaries that are pseudo-10-locular and, therefore, resemble the 10-locular ovaries of 
Gaylussacia (Vander Kloet & Dickinson 1992). 

<a>-' 

Gaylussacia usually has been divided into three sections (following Sleumer 1967): section Vifis-uliiflj 
Drude (only G. brachyccra), section Gaylussacia (ca. 47 species, mainly South American; plants mostly evei 
green, with stalked glands on the leaves), and the North American section Dccameriutn Torr. & A. Gray (a 
live species; plants deciduous, with sessile glands on the leaves). The last is the focus of this study. All three 
sections are represented in eastern North America: G. brachyccra, four of the species of section Guy/ussaciu 
(G. mosicn Small, (.. ciumosa (Andr.) A. Gray, G. orocola (Small) Camp, G. bigeloviana (Fernald) Sorrief 




Weakley; see Sorne & Weakley, 2007, although the latter two entities are often included within an exp< 

G. dumosa), and all five species of section Dccamcrium [G. frondosa (L.) Torr. <Sr A. Gray, G. tomcnum (A. 
Gray) Pursh ex Small, G. nana (A. Gray) Small, G. ursinu (Curtis) Torr. &r A. Gray, and G. baccaia (Wang.' K 
Koch], For more information on the history and systematics of these groups see Camp (1935), Wood (1961). 
Luteyn et al. (1996), Floyd (2002), and Sorrie and Weakley (2007). 

A phylogenetic analysis of Gaylussacia had not been carried out prior to that of Floyd (2002). In her 
morphology-based analysis, which included most of the species of the genus, the monophyly of sects. Deuimt- 
rium and (of course) VUis-idaea were supported, but sect. Gaylussacia was paraphyletic. All three traditional 
sectional divisions were weakly to moderately supported by cpDNA (ti nL-trnF) data, as well as in analyses 
combining morphological and molecular data, but nrlTS sequences alone did not support the mono 
of either sect Decamcrium or sect. Gaylussacia. Molecular data were lacking for most species of sect. Gaylus¬ 
sacia and relationships within this large clade remain poorly understood. Floyd’s (2002) analysis also did 
not clarify the placement of G. brachyccra', neither did it fully clarify species delimitations and relationship 
within sect. Decamcrium. 

Camp (1941) recognized three subsections within sect. Decamcrium, two ol which, subsects. Baud.* 
and Ursinac. are monotypic, i.e., Gaylussacia baccaia and G. ursinu, respectively. The remaining three W*> 
i.e., G. fnvidasa, G. lomentosa and G. nana, were placed in his subsection Frondosae. The taxa within the 
Frondosac group have had a confusing taxonomic history and have been variously circumscribed. Ellie 
(1821) and Chapman (1889) recognized only G. frondosa. Radford et al. (1964) recognized only G J J 
but included two varieties: i.e., G. frondosa var. tomento.su A. Gray and var .Jrondosa\ Wunderlin and Hansen 
(2003) did the same, recognizing only var. tomentosa as occurring in Florida. Gray (1878) was the first4® 
thor to recognize G. Jrondosa var. tomontosu, and he also recognized G. frondosa var. nana A. Gray Harpe' 
(1906) treated the members of subsect. Frondosae as a single species with three varieties, i.e.. G. 
var frondosa, var tomerUosa, and \ar. nana. Sleumer (1967) also treated the subsection as a single specif" 
but recognized O. ffondosa var. polycodioidcs Camp and f. glaucophylla Camp (both pertaining to norths® 
plants usual 1> treated within var. frondosa). Small (1897, 1933), Camp (1935, 1941), Wood (1961). Dunctf 
and Brittain (1066), and Luteyn et al. (1996) treated the subsection as three separate species (G.j condor- 
tomentosa. and G. nana), Most recently, Floyd (2002) treated the subsection as one species with three whet* 
ies. However, she suggested that var. nana perhaps should be considered a separate species, citing a possM 
lineav orting oi hybridization event as being responsible for its incongruent placement in her nrlTS * 
cpDNA trees. \et, as noted by Floyd, such problems could also merely be due to inadequate data, 
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resulting in poor resolution in her eladograms. She recommended more research in order to resolve species 
limits within the Frondosae group. 

Field observations (over many years, by W.S. judd) of the notable and seemingly consistent morphologi¬ 
cal differences between the sympatric Gaylussacia tomcntosa and G. nana in Florida motivated the present 
study (Fig. 1) The purpose of this study is to resolve the phylogenetic relationships within this subsection, 
i.e., the G. frondosa complex, by focusing on plants identified as G. frondosa, G. tomento$<t and G. nana. We 
sought to determine appropriate species limits within the G. fronftosd complex, based on phenetic, evolu¬ 
tionary* diagnostic, and apomorphic species concepts (Davis Sr Nixon 1992; Donoghue 1985; de Queiroz 
07; Judd et al. 2007; Mishler 1985; Mishler & Theriot 2000; Wheeler Sr Platnick 2000; Wiley & Mayden 
2000 ). 


MATERIALS & METHODS 

Taxon sampling and field work, —Voucher material representing the Gaylussacia frondosa complex was col¬ 
lected by M.T. Gajdeczka and W.S. Judd from numerous localities in Florida, Georgia, South Carolina, and 
orth Carolina in the spring and summer of 2007 in order to estimate genetic diversity (Table 1; all deposited 
at FLAS). Gaylussacia Ursirui and G. bacccita were also collected, and were selected as outgroups based on 
ie traditional classification of section Occam tfkitn (Camp 1941) and the recent phylogeny of Floyd (2002). 
Leaf material for DNA extraction of field-collected specimens was preserved in silica gel. DNAs derived 
from herbarium collections also were included in the study to supplement the number of evolutionary units 
in the analyses, especially of G. f rondosa , as well as outgroup taxa. One collection of G. dumosd (of section 
Gaylassada) was included as a more distant outgroup. 

Morphology. —Potentially phylogenetically informative morphological characters were selected after 
‘.areful consideration of the pattern of variation seen in material at the University of Florida Herbarium of 
die Florida Museum of Natural History (FLAS)* field observations (of M .T Gajdeczka and W.S. Judd), and 
previous taxonomic work on the genus (Duncan &r Brittain 1966; Luteyn et al. 1996; Floyd 2002). Seventy 
qualitative and quantitative morphological characters were initially assessed; all above ground organs were 
examined, including both vegetative and reproductive features, with an emphasis on density, distribution, 
J nd form of the hairs (both glandular and non-glandular). From the initial list of characters observed, 20 
were selected as potentially phylogenetically informative and were included in the matrix for analysis (see 

Tables 2,3). 

Unicellular hair presence, density and length and resin gland presence, density and width have been 
'nsidered especially useful in distinguishing species of Gqylussarici (Floyd 2002; Luteyn et al. 1996), and 
ui$ these were included in the analyses. Most of the initial hair and gland measurements could be grouped 
lf uo character-clusters, which showed identical or very similar patterns of variation among taxa. In such 
a single character (from each group of highly correlated characters) was selected for inclusion in the 


bus 



^netic analyses. However, measurement of the length of the longest hair and average hair length of 
e adaxial leaf surfaces (chars. #5 and 6; see Table 2) were both included in order to take into account the 
r \ed variation in the range of hair length between specimens. 

Many morphological characters used in the analyses (Table 2) were readily divisible into discrete states, 
aiding arbitrary decisions relating to state delimitation (Stevens 1991V However quantitative characters 
more problematic. Variation in these characters (e.g., chars. #1, 3, 6, 7,8. 11, 12. and 13, see Table 2) 
messed by means of bar graphs, and the states of those included in the analyses were delimited by more 
* ess discrete gaps (e.g., Fig. 2, char. #9). The most problematic characters in terms of state delimitation 
^an-glandular hair length on the abaxial leaf surface (#8), and the length (#33) and width (#34) of the 
gest blad es; these characters showed nearly continuous distributions, but were included in the analyses 
*R USe ^ ave been stressed as taxonomically important in the Q<iyluss(ici&.frofidosA complex (Duncan 
ittain 1966), and because the overlap in values between taxa was limited. Many characters could not 
luded in the analyses because they showed too much infraspecific variation or varied continuously 
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Gaylussacia Kunih is a New World genus of 53 species of mainly undersiory shrubs occurring in mesicto 
xeric woodlands and shrublands, as well as acidic bogs. The genus is differentiated from the closely-related 
Vaccinium L. by ten-locular ovaries, drupaceous fruits containing ten pits, the presence of resin glands on the 
leaves (except for G. brachycera (Michx.) A. Gray), and the lack of staminal spurs (Duncan &r Brittain 196' 
Luteyn et al. 1996; Palser 1961; Wood 1961). However, the phylogenetic relationships between members of 
the two genera are still uncertain (Kron et al. 2002b). Vacarmm was shown to be non-monophyletic as tradi¬ 
tionally circumscribed, with Gaylussacia nested within it (Kron et al. 2002a). These results are corroborated 
by the ambiguous generic identity of G. brachycera , which, in addition to being eglandular and similar in 
form to Vaccinium species, has been shown to be sister to the rest of Gaylussacia (Floyd 2002). Finally, some 


Vaccinium species have ovaries that are pseudo-10-locular and, therefore, resemble the 10-locular ovariesol 


Gaylussacia (Vander Kloet <$r Dickinson 1992). 

■Jr 

Gay/uss ,acia usually has been divided into three sections (following Sleumer 1967): section Vitis-iiaw 
Drude (only G. brachycera'), section Gaylussacia (ca. 47 species, mainly South American; plants mostly eur- 
green, with stalked glands on the leaves), and the North American section Dccamerium Torr. & A. Gray tea 
five species; plants deciduous, with sessile glands on the leaves). The last is the focus of this study. All three 
sections are represented in eastern North America: G. brachycera, four of the species of section Qjyfusstuiu 
iG. mosieri Small, G. dumosa (Andt.) A. Gray, G. orocola (Small) Camp, G. bigehniana (Fernald) Sorrie& 
Weakley; see some &r Weakley, 2007, although the latter two entities are often included within an expanded 
G. dumosa), and all hve species of section Dccamerium JG. frondosa (L.) Torr. & A. Gray, G. tomentosa (A 
Gray) Pursh ex Small, G. nana (A. Gray) Small, G. ursina (Curtis) Torr. & A. Gray, and G. baccata (Wang.)K. 
Koch). For more information on the history and systematics of these groups see Camp (1935), Wood (1961 


Luteyn et al. 11996), Floyd (2002), and Sorrie and Weakley (2007). 

A phylogenetic analysis of Gaylussacia had not been carried out prior to that of Floyd (2002). In her 
morphology-based analysis, which included most of the species of the genus, the monophyly of sects. Deaw i 
rium and (of course) Vitis-ulucii were supported, but sect. Gaylussacia was paraphyletic. All three tradition.!, 
sectional divisions were weakly to moderately supported by cpDNA ( tmL-trnF ) data, as well as in analyse? 
combining morphological and molecular data, but nrlTS sequences alone did not support the mono 
ol either sect. Decamerium or sect. Gtitylussttcid. Molecular data were lacking for most species of sect. Gaylu 
sacia and relationships within this large clade remain poorly understood. Floyd's (2002) analysis also did 
not clarify the placement of G. brachycera', neither did it fully clarify species delimitations and relationship 
within sect. Decamerium. 

Camp i 1941) recognized three subsections within sect. Decamerium, two of which, subsects. Bat :e 






and L/rsinae, are monotypic, i.e., Gaylussacia baccata and G. ursina, respectively. The remaining three tax* 
i.e., G. frondosa. G. lomentosa and G. nana, were placed in his subsection Frondosae. The taxa within k. 
Frondosae group have had a confusing taxonomic history and have been variously circumscribed. Elliott 
1 1S 21 ) and Chapman (1889) recognized only G. frondosa. Radford et al. (1964) recognized only G. frontlet 
but included two varieties: i.e., G. frondosa var. tomentosa A. Gray and var. frondosa; Wunderlin and Han* 
2003) did the same, recognizing only var. tomentosa as occurring in Florida. Gray (1878) was the first« 
thor to recognize G. frondosa var. tomentosa, and he also recognized G. frondosa var. nana A. Gray. Harpef 
(19061 treated the members of subsect. Frondosae as a single species with three varieties, i.e., G.pon 
var. frondosa, var lomentosa, and var. nana. Sleumer (1967) also treated the subsection as a single specit 
but recognized G. frondosa var. polycodioides Camp and f, glaucophvlla Camp (both pertaining to northed 
plants usually treated within var. frondosa). Small (1897, 1933), Camp (1935, 1941), Wood (196P. Dune; 
nd Brittain (1966'. and Luteyn el al 11996) treated the subsection as three separate species (G.frtMaOsa, 
ti>mcntosa. and G. nana). Most recently. Floyd (2002) treated the subsection as one species with three var' 

* * r 

it‘s However, she suggested that var. nana perhaps should be considered a separate species, citing a 
lineage sorting or hybridization event as being responsible for its incongruent placement in her nr!T. s 1 
cpDNA trees. Yet, as noted by Flo\J, such problems could also merely be due to inadequate data. 
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resulting in poor resoluiion in he r cfodograms. She recommended more research in order 10 resolve species 

s within the Frondosae group. 

Field observations (over many years, by W . Judd) of the notable and seemingly consistent morphologi¬ 
cal differences between the sympatric Gavlusscuia lomcntosa and G. nurui in Florida motivated the present 
study (Fig. 1). The purpose of this study is to resolve the phylogenetic relationships within this subsection, 
U.,the G. frondosa complex, by focusing on plants identified as G. frondosa, G. tomentosa and G. nana. We 
sought to determine appropriate species limits within the G. frondosa complex, based on phenetic, evolu- 
onary, diagnostic, and apomorphic species concepts (Davis & Nixon 1992; Donoghue 1985; de Queiroz 
AW; Judd et al. 2007; Mishler 1985; Mishler & Theriot 2000; Wheeler & Platnick 2000; Wiley & Mayden 
2000 ). 


MATERIALS & METHODS 


Taxon sampling and field work .—Voucher material representing the Gaylussada frondosa complex was col- 
dbyM.T. Gajdeczka and W.S. Judd from numerous localities in Florida, Georgia, South Carolina, and 
North Carolina in the spring and summer of 2007 in order to estimate genetic diversity (Table 1; all deposited 
at FLAS). Goyfussacia ursina and G. baccala were also collected, and were selected as outgroups based on 
'he traditional classification of section Decameriutn (Camp 1941) and the recent phylogeny of Floyd (2002). 
leaf material for DNA extraction of field-collected specimens was preserved in silica gel. DNAs derived 
from herbarium collections also were included in the study to supplement the number ol evolutionary units 
in the analyses, especially of G. frondosa, as well as outgroup taxa. One collection of G. dumosa (of section 



was included as a more distant outgroup. 

Morphology .—Potentially phvlogeneiically informative morphological characters were selected after 

ireful consideration of the pattern of variation seen in material at the University of Florida Herbarium ol 

ie Florida Museum of Natural History (FLAS), field observations (of M.T. Gajdeczka and W.S. Judd), and 

previous taxonomic work on the genus (Duncan & Brittain 1966; Luteyn et al. 1996; Floyd 2002). Seventy 

qualitative and quantitative morphological characters were initially assessed; all above ground organs were 

examined, including both vegetative and reproductive features, with an emphasis on density 7 , distribution, 

md form of the hairs (both glandular and non-glandular). From the initial list of characters observed, 20 

were selected as potentially phylogenetically informative and were included in the matrix for analysis (see 
tables 2,3). 

Unicellular hair presence, density and length and resin gland presence, density and width have been 

considered especially useful in distinguishing species ot Gttyhlssactil (Floyd 2002; Luteyn et al. 1996), and 

u s these were included in the analyses. Most of the initial hair and gland measurements could be grouped 
k c sracter-clusters, which showed identical or very similar patterns of variation among taxa. In such 

# J I 

s - a single character (from each group of highly correlated characters) was selected for inclusion in the 


Phylogenetic 



analyses. However, measurement of the length ol the longest hair and average hair length of 
surfaces (chars. #5 and 6; see Table 2) were both included in order to take into account the 
' ed var lation in the range ot hair length between specimens. 

y tnorphological characters used in the analyses (Table 2) were readily divisible into discrete states, 
^ding arbitrary decisions relating to state delimitation (Stevens 1991). However quantitative characters 
r edmore problematic. Variation in these characters (e.g., chars. #1, 3, 6, 7, 8,11, 12, and 13; see Table 2) 
0r | 3SSesscd fe y m eans of bar graphs, and the states of those included in the analyses were delimited by more 
discrete gaps (e.g., Fig. 2, char. #9). The most problematic characters in terms of state delimitation 
^ g andular hair length on the abaxial leaf surface (#8), and the length (#33) and width (#34) of the 

ades, these characters showed nearly continuous distributions, but were included in the analyses 

uccdusg t u ^ i * 

^ have been stressed as taxonomically important in the Guylusscictd frondosa complex (Duncan 
lain 1966), and because the overlap in values between taxa was limited. Many characters could not 

c inchiHori % r 

tn the analyses because they showed too much infraspecific variation or varied continuously 
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Fig. 1. Gaylussacia nana (top of A, B ( D) and 6. tomentosa (bottom of A r C, E). A. The glaucousness (gray coloration) of the abaxial leaf surface of G. m - 
compared to G. tomentosa. B&C. SEM image of glaucous abaxial surface of the leaf of 6. nana (B) and non-glaucous abaxial leaf surface of 6. foment 
(C). D. Habit of G. nana , E. habit of 6. tomentosa , Photographs A, D and E, from Morningside Nature Center, Alachua Co., FI, by W.S. Judd; photograp 
BandCby G.M. lonta. 


s wert 


across taxa (and thus could not be delimited into states). Characters not observed for particular species 
scored as missing values, as were situations where a character was considered “not applicable” in a particular 
taxon, and coded as u ? w (table 2). 

DNA extraction , amplification and sequencing . — Leaves and flowers dried in silica-gel were used for 
DNA extraction of field collected material. A modified version of the 2x CTAB procedure (Doyle & 
1987) was used, as described in Whitten el al. (2007) with the addition of proteinase K (5 units) to t 
extraction buffer instead of 2-mercaptoethanol. Precipitated DNA pellets from extractions of field-collected 
material were resuspended lit 200 pL of Tris-EDTA (TE) buffer. Small fragments of leaf tissue were remo' 
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iutSpecimens sampled for the morphological and nrlTS analyses. All vouchers deposited at FLAS. 
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Voucher 

Collection Location 

ITS-GenBank 

4. l/oscoff 405 

King & Queen Co., VA 

NA 

EtV Wood 4222 

Clinton Co., PA 

FJ985237 

W.S. Judd 1657 

Barnstable Co., MA 

NA 

W.S. Judd 716 

Norfolk Co., MA 

FJ985236 

M. Gajdeczka 50 

Burke Co., NC 

FJ985196 

M. Gajdeczka 18 

Alachua Co., FL 

FJ985199 

A.B. Pittman 05090108 

Orangeburg Co., SC 

FJ985228 

AS. Pittman 053009303 

Aiken Co., SC 

FJ985227 

G.W. Seckinger, Jr. 271 

Georgetown Co., SC 

FJ985226 

J.B. Nelson 16473 

Richland Co., SC 

NA 

M. Gajdeczka 56 

Colleton Co., SC 

FJ985197 

M. Gajdeczka 57 

Colleton Co., SC 

FJ985212 

R.R. Smith 1563 

Santa Rosa Co., FL 

NA 

W.S. Judd 786 

Norfolk Co., MA 

NA 

C. van Hoek & B. Wargo 1084 

Polk Co., FL 

FJ985225 

J. Amoroso 261 

Levy Co., FL 

FJ985223 

M. Gajdeczka 16 

Alachua Co., FL 

FJ985191 

M. Gajdeczka 22 

Putnam Co., FL 

FJ985202 

M. Gajdeczka 23 

Suwannee Co., FL 

FJ985203 

M. Gajdeczka 31 

Lowndes Co., GA 

FJ985206 

M. Gajdeczka 38 

Jeff Davis Co., GA 

FJ985194 

M. Gajdeczka 60 

Alachua Co., FL 

FJ985213 

M. Gajdeczka 61 

Alachua Co., FL 

FJ985214 

M. Gajdeczka 64 

Alachua Co., FL 

FJ985217 

M. Gajdeczka 65 

Alachua Co., FL 

FJ985218 

M. Gajdeczka 66 
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Longest Abaxial Leaf Hair Length 


G. baccata 


G. ursina 


G. frond os a 


G. nana 


G. tomentosa 



Average Length (mm) 


Fig. 2 . Graph showing variation in the length of the longest non-glandular hair on the abaxial surface of the leaf for members of Gaylussacia sect 

Decamerium. Squares represent the average of the measurements for all specimens of the corresponding taxon. Lines on either side represent the 
range of all measurements. 


1m 2. Morphological characters and character states included in the phylogenetic analyses of Gaylussacia. 



1. 

2 
A— * 

5. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10. 
11 . 

12 . 

13 . 

14 

15. 

16. 

17 . 

18 . 

19 . 

20 . 


Plant height: n < 0.7 m (0); n > 0.7 m (1). 

Ratio of average length of lower secondary branches to length of main stem: n < 7 (0); n > 7 (1). 

Length of longest non-glandular hair on branch: n < 0.5 mm (0); n > 0.5 mm (1). 

Number (density) of non-glandular hairs on the adaxial leaf surface per standardized 2.56 mm 3 region: none (0); 0 < n 

< 10(l);n> 10(2). 

Length of longest non-glandular hair on adaxial leaf surface: 0 < n < 0.25 mm (0); n > 0.25 mm (1). 

Length of average non-glandular hair on adaxial leaf surface; 0 < n < 0.175 mm (0); 0.175 < n < 0.322 mm (1); n 2 0322 
mm (2). 

jm 

Number (density) of non-glandular hairs on the abaxial leaf surface per standardized 1.32 mm 3 region: 0 < n < 6.5 (0 
<n<40(l);n>40(2). 

Length of longest non-glandular hair on abaxial leaf surface: 0 < n < 0.275 mm (0); 0.275 < n < 0.45 mm (1); n > 0.45 mm (2' 

Number (density) of resin glands on abaxial leaf surface per standardized 4.0 mm 3 region: 0 < n < 19 (0); 19 < n < 65 (l 
n > 65 (2). 

Glaucousness of leaves: none (0); slight (1); moderate (2); strong (3). 

Length of blade (average of five longest leaves per specimen): n < 41.1 mm (0); 41.1 < n < 70 mm (1); n > 70 mm (2). 
Width of blade (average of five longest leaves): n < 20.3 mm (0); 20.3 < n < 30 mm (1); n > 30 (2). 

Petiole length (average of five longest leaves): n < 1.35 mm (0); n > 1.35 mm (1). 

Bract type: leaf-like (0); small (1; states coded and scored as per Luteyn et al. 1996, and for this character and #16, 1 
and 19, our observations, when these were possible, always matched these scorings). 

Resin glands on adaxial leaf surface: present (0); absent (1). 

Fruit color: black (0): blue (1; states coded and scored as per Luteyn et al. 1996). 

Corolla color: white (0); greenish-white to pinkish white (1); red (2; states coded and scored as per Luteyn et al. 1996 , 1 
Corolla shape: cylindrically urceolate (0); broadly urceolate (1; states coded and scored as per Luteyn et al 1996). 

Pilose hairs on the filaments: present (0); absent (1; states coded and scored as per Luteyn et al. 1996). 

Exudate around perimeter of resin glands on leaves: present (0); absent (1). 
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specimens tor DNA extraction. Anticipating a lower DNA concentration, the herbarium 

m were ^suspended in only 120 pL of Tris-EDTA (TE) buffer. In all extractions the resuspended 

purified using Qiaquick columns (Oiagen Valencia, California, USA) and Buffer PE, then eluted 
ulfer EB. 

s were carried out using Sigma Jumpstart Tar/ polymerase and reagents (Sigma-Aldrich, 
i EISA). Initially a survey of phylogenetically informative regions were carried out using 
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S region and seven plastid DNA regions ( trnL-F, rpl32-trnL, trnQ-rpsl6 , tmH-psbA , 
PsbD-trnT and trnV-ndhC) based on suggestions of phylogenetically useful regions in Small et al. 
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(1998) and Shaw et al. (2005, 2007). Two representative field-collected specimens of each of the three in¬ 
group taxa ( Gaylussacia frondosa, G. tomentosa, and G. nana) and one sample each of G. baccata and G. ursina 
were included in the survey. The reaction mixture tor amplification of nrlTS from field collected specimens 
included 7.0 pL of betaine (5 M), 12 pi. 11:0, 2.5 pL 10X buffer, 2.0 pL MgCh (25 mM), 0.5 pL dNTPsUO 
pM), 0.5 pL of each primer (10 pM), 0.5 pL template, and 0.2 pL Taq polymerase (25.5 pL total). In order 
to improve amplicon concentration, volumes were adjusted to 1.0 pL template and 11.0 pL H 2 O (25.7 pL 
total). Reaction mixtures for the amplifications from herbarium specimens totaled 24.7 pL (18 pL H.O,2.5 
pL 10X buffer, 2.0 pL MgCh (25 mM), 0.5 pL dNTPs (10 pM), 0.5 pL each of forward and reverse primer 
(10 pM), 0.5 pL template, and 0.2 pL Taq polymerase). 

The following PCR protocol was used for the initial survey: an initial bake at 94°C for 3 min; then 30 
cycles of (1) denaturation at 94°C for 30 s, (2) annealing at 60°C for 30 s, and (3) extension at 72°C for 2 
min; and a final extension at 72°C for 3 min. Primers 17SE and 26SE from Sun et al. (1994) and the PCR 
program in Whitten, et al. (2007) were used for amplification and sequencing of the nrlTS region from field 
collected specimens. With the aim of improving nrlTS amplicon yield from herbarium specimens, primer- 
ITS A/ITS C and ITS B/ITS D and the protocol from Blattner (1999) were used, and the reaction mixture was 
doubled to 50.4 pL (30 pL HiO, 9 pL MgCh (25 mM), 5 pL 10X buffer, 4 pL template 2 pL each of forward 
and reverse primer (lOpM), 1.0 pL dNTPs (lOpM), and 0.4 pL Taq polymerase). 

Phylogenetic analyses .—Herbarium specimens, based on the first author's own field work or Irom 
the collections of FLAS, and representing populations of all the entities within the Gaylussacia frondosa 
complex, were employed as terminal taxa in the phylogenetic analyses. These analyses served as the basb 
not only of the assessment of phylogenetic relationships within Gaylussacia section Decamerium but also 
for hypotheses of species circumscription (based on an application of the phylogenetic/apomorphic speck 
concept; Donoghue 1985; Mishler 1985; Mishler & Theirot 20001. 

PAUP* 4.0b 10 tSwofford 2002) was used to construct most parsimonious trees for three data sets 1) 
a morphological analysis of 38 specimens (i.e., terminal taxa, representing populations; 26 of these field- 
collected and 12 herbarium/FLAS collections) and 20 characters, and 19 of these parsimony informative 

0 - A 

(see Table 2), (2) a nrlTS analysis of 47 specimens and 827-857bp (31 field-collected specimens) or 648-otrt 
bp (FLAS herbarium material, 16 specimens) but only 12 parsimony-informative sites (and 19 total variable 


characters), and (3) a combined "total evidence" analysis, i.e., morphology + nrlTS, for a pruned data set of 
31 specimens (24 field-collected, and 7 herbarium), and 48 variable characters, with 31 of these parsimony 
informative. All characters were equally weighted and unordered. The morphological analysis used a maxi¬ 
mum parsimony (MP) heuristic approach, 500 random-addition replicates, TBR, MaxTree = 200 per replicate, 
and MulTrees on. Relative support for clades in all analyses was evaluated using a fast-heuristic bootstrap 
analysis (1000 replicates in the morphological analysis, 100 replicates in the nrlTS and combined analyst 
The nrlTS and analysis used a MP heuristic approach, 100 random-addition replicates, TBR, MaxTree = 
replicate, and MulTrees on. The combined analysis was similar but used MaxTree = 10,000. Gaps were codeu 
as missi ng data; indels were not coded as characters. Morphological character state changes were mapped onto 
the morphological cladogram and traced onto the total evidence MP cladogram using the PAUP* output» 
MacClade 4.05 (Maddison &r Maddison, 2005). The molecular and morphological analyses resulted in trees 
with no strongly supported incongruent patterns, so were combined (see “total evidence" analysis, abov. 



RESULTS 


Morphological analysis .—A heuristic search using maximum parsimony yielded 3360 most parsmiomo’ 
trees (MPTs) u ith a length of 65 (Consistency Index |CI] = 0.462, Retention Index [Rll = 0.829, Rescaled 
Consistency Index (RC] = 0.383). The strict consensus tree is shown in Figure 3. Gavlussncia tomentosa. v 
nana. G. frondosa , and G. lirsinu form a clade with moderate support (bootstrap [BS] = 78%) in all trees (Fi^ 


3A) The G. frondosa complex, i.e., G. tomentosa, G. nana, and G. frondosa, also forms a clade (BS = 85%: *1 
3). Within the G. frondosa complex, populations of G. tomentosa form a clade (without bootstrap supp orl 
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A) Morphology 
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G, tomentosa MG30 
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G. tomentosa MG36 
G. tomentosa MG3 7 
G. Comemosa MG5S 
G. tomentosa x nana MG43 
G. ursina AC526* 

G. ursina MG47 
G. ursina MG48 
G. ursina AD 1909 
G. baccata MG50 
G. baccata \Vf7l6 
G. baccata EW4222 
G. baccata AV405* 

G. baccata WJ1657* 

G. dumosa MG 18 


8) nrlTS 
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-G. tomentosa $022613* 

- G. tomentosa BT698* 

--- G. tomeniosaJN2Q50Q* 

■ G. frondosa GS271 

-GJroruiosa AP05309303* 

- G f i imdosa AP05090 108 

-—- G. frondosa MG56 

‘-—-— G. frondosa MG57 

- G. nana MGJ 6 

— --—-— G. nana MG38 

--—- G. nana MG22 

-—-G. nano MG23 

---- G. nana MGJi 

-*-G. nana MG60 

- G. nana MG6I 

I--—- G. nana MG64* 
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-G. nana MG66" 

-— G. nana MG70* 
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-- G. nana]A261 * 

1--— -G. nana CVH1 084 * 
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- G baccata MG50 
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ht 


^ ^ con sensus trees resulting from (A) the morphological analysis and (B) the nrlTS analysis. Collector's names abbreviated by initials (see Table 

jWfllfliHr analysis indicated with an asterisk. See results for bootstrap values. 


sister to a clade comprising the populations of G. frondosa and G. nana (without bootstrap support) 


u jrondosa + G. nana clade. the populations of G. nana form a clade (BS = 58%), which is nested 
Paraphyletic G. Jrondosa. The terminals representing populations of G. urshiu also form a clade (BS 
vaccata is unresolved at the base of the Decamerium clade. 
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nrlTS analysis.—None of the plastid regions showed phylogeneiically intormalive variation, so after 
preliminary work these regions were not pursued. However, the 47 sequences of the nrlTS region did show 
useful variation and, therefore, were studied further. A MP heuristic search (100 replicates, MaxTrees = 10 
per replicated of the nrlTS sequences resulted in a total of 1000 MPTs with a length 33 (Cl = 0.970, RI = 
0.992, RC = 0.962). The strict consensus tree is shown in Figure 3B. The monophyly of Gaylussacia nana 
(BS = 56%) was weakly supported. A G. tomcntosa + G. frondosa clade (BS = 95%) was strongly supported 
A third clade in the strict consensus tree groups G. baccata and G. ursina (BS = 69%). 

Combined analysis. —The combined total evidence analysis included all specimens that had both 
successfully sequenced nrlTS regions and morphological data. A MP heuristic search resulted in a total ol 
1000 MPTs of length 89 (Cl = 0.663, RI = 0.878, RC = 0.582). A representative tree is shown in Figure 4, 
with the strict consensus and bootstrap consensus values mapped. Clades Gaylussacia nana (BS = 91%) and 
G. tomcntosa + G. frondosa t BS = 87%) are strongly supported and are sister in the strict consensus nee, thu 
supporting the monophyly of the G. Jrondosa complex (BS = 62%). All specimens of G. tomcntosa lorm t 
monophyletic group in the strict consensus tree (but without BS support). Specimen Gajdtxzka 43 (MG4 
probably a hybrid (see below), and the operational taxonomic units (OTUs) representing G. frondosa forma 
paraphyletic grade in the G. tomcntosa + G. frondosa clade, and G. frondosa forms a clade only in some of the 
equally parsimonious trees. The G. tomcntosa + G. frondosa + G. nana clade is sister to a strongly supported 
G. ursina clade in the SC (BS = 95%). A G. baccata clade (BS = 100%) is sister to the G. ursina + G. tomcntosa 
+ G. Jrondosa + G. nana clade (BS = 59 

Morphological synapomorphies. —Mapping morphological characters onto a randomly chosen MP I 
from the combined (or total evidence) analysis showed that most clades were supported by at least two mor¬ 
phological synapomorphies. The monophyly of Gaylussacia nana in the combined MPT was supported by six 
character states: short non-glandular hairs on both adaxial and abaxial surfaces of the leaf (chars. #5,6. and 
, Table 2), short and narrow blades (#11 and 12) and short petioles (#13). The monophyly of G. tofnentosa uas 
supported by dense non-glandular hairs on both adaxial and abaxial leaf surfaces (chars. #4 and 7), as well 
as long hairs on the branches #3) and on the abaxial leal stir I ace (#8'. Gavlussm id tomcntosa specimen sa»o 
were less glaucous than those of G. frondosa and especially G. nana (char. #10, see also Fig. 1). Gaylussac 
Jrondosa (when monophyletic) was poorly' supported by morphological synapomorphies: most G. frondosa 
specimens are characterized by sparse or absent hairs on the adaxial leal surface (char. #4) and mediun 
length leaf hairs on the abaxial surface (char. #8), although the variation pattern is homoplasious. In addition 
Gaylussacia frondosa tends to have fewer unicellular hairs on its twigs than either G. tomcntosa or G. tw° d 

Specimens of Gaylussacia ursina formed a monophyletic group supported by the average length of the 
non-glandular hairs on the adaxial leaf surface (char. #6) (as opposed to the 'longest hair” measurement, 
long abaxial leaf hairs (#8), the longest and widest blades of all taxa in the analysis (#11 and 12), usual!) 
nonglaucous It ives (#10), and sparse glands on the abaxial leaf surface (#9). The monophyly of G. bate 
was supported by five unique synapomorphies: very dense abaxial leaf glands (char, #9), a red corolla (#” - 
a cylindrically urceolate corolla (#18), and presence of glands on the adaxial leaf surface (#1 n Addition! 
derived features that distinguished it were sparse adaxial leaf hairs (char. #7), medium length abaxuii tea 

hairs (#8), black coloration of truits (#16), presence of filament hairs (#19), lack of glaucousness (#10), an 
narrow leaf blades (#12). 

The GaWussaria tomcntosa + G. frondosa clade was supported mainly by features that differentiate it 
from G. nanu: tall plant height (char #1), long adaxial leaf hairs (#5 and 6), dense abaxial leaf glands (# 
long and wide leaf blades (#11 and 12), and long petioles (#13). The monophyly of the G. tomentosu + d 
jrondosa + G. nuna clade was supported by dense abaxial leaf hairs (char. #30), leaves with slight to intend 
glaucou>ness (#10, in many specimens, but lost in some, especially G. tomcntosa), blue fruits (#16), andgb 
brous filaments (#19). Finally, the monophyly of the G. tomentosa + G. frotulosa + G. nana + G. u/sina 


f 

t 


was supported by a lack of adaxial leaf glands (char. # 15 ), greenish-white to pinkish-white coloration v 
of the broadly urceolate corolla (#18). and resin glands lacking exudate (#20). 
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G. dumosa MG 18 
G. tomentosa WJ1679 
G. tomentosa MG32 

G. tomentosa MG33 
G. tomentosa MG37 

G. tomentosa MG20 

G, tomentosa MG58 
G. tomentosa MG24 
G. tomentosa MG30 
G. tomentosa MG17 
G. tomentosa MG36 
G. tomentosa MG 19 
G. frondosa AP05090108 

G. frondosa MG57 

G. frondosa GS271 

G. frondosa MG56 

G. tomentosa x nana MG43 

G. nana MG 16 

G. nana MG22 

G. nana MG23 

G. nana MG61 

G. nana RG8484583 

G. nana MG38 

G. nana MG60 

G. nana MG31 

G. ursina ADI909 

G. ursina MG47 
G. ursina MG48 



G. baccata MG50 
G. baccata WJ716 
G. baccata EW4222 


tree from the combined nrlTS and morphology analysis. The numbers above the lines represent the number of character changes. 
^abbreT" 1 ^ Stfitt consensus tree are indicated by bold lines; bootstrap consensus values provided below the lines. Individual OTUs indicated 
ation of collector names along with the collection number (see Table 1). 
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Putative hybrids .—V arious specimens are putatively hybrids between G. tomentosa and G. nana f'O ryfl 
& Bridges 19346, Herring 1548, Gajdeczka 21, Gajdeczka 43, Gajdeczka 75, and Gajdeczka 76) or between G. 
frondosa and G. nana (Judd 3118) and all are deposited at FLAS. Most of these specimens were excluded from 
analysis because of the intermediate nature of the data. The interspecific hybrid status of these specimens is 
supported by both DNA and morphology. The putative G. tomentosa x nano hybrids are noteworthy in having 
leaves that are usually glaucous (and thus similar to G. nano) but also have long unicellular hairs on their 
abaxial surfaces (like G. tomentosa ); their leaves often are also smaller then is characteristic of G. totnentoso 
(especially so in Gajdeczka 21), again like G. nana. The specimen Judd 3118 has small, abaxially glaucous 
leaves (similar to those of G. nana), but their adaxial surfaces are essentially glabrous (like G. frondosa), and 
the twigs are only sparsely pubescent (again, similar to G. frondosa). For nrlTS, there exist polymorphisms 
at all of the variable nucleotide positions among these three species. For example, at base 163 in our ITS 
matrix, there is a variable nucleotide C or T, with C occurring in G. tomcntosa/jrondosa vs. T in G. mim 



hybrids exhibit a strong poly morphic signal for both C and T. Additional polymorphic sites (substitution' 


arc found at sites 207, 344, 358,679, 682, as well as two indels (insertion/deletion) from bases 704-711 



a single base indel at 739. Because we did not clone these PCR products to separate the different copies ol 
ITS, we cannot be certain of the distribution of these characters and their association with either parental 
species. However, it has been shown that PCR-mediated recombination can cause false or mixed signal of 
the original parental types (Lahr &r Katz 2009). Nonetheless, we are confident that these sequences represent 
interspecific hybrids because of the high association of the polymorphic nucleotide positions only at the 
variable sites among different species and at every variable site among the species in question (see Soltis el 
al. 2008) and because of the intermediate morphology discussed above. 


in 


DISCUSSION 

Taking all three analyses into account, it is clear that Gaylussacia nana and G. tomentosa are eladospecies 
(Donoghue 1985; Mishler 1985; Mishler &r Theriot 2000), i.e., each of these species is a monophyletic group 
of populations (see Donoghue 1985; Mishler 1985; Mishler &r Theriot 2000). Thus, both should be recognized 
as distinct species. We note that they are largely sympatric and are easily distinguished by the morphologi¬ 
cal characters listed above. The monophyly of G. nana is supported by all analyses, with strong support 
the combined analysis (see Fig. 4, BS = 91%). Diagnosable features for G. nana included short hairs on both 
adaxial and abaxial leaf surfaces, short and narrow blades, and short petioles (and compare these with 
features of G. tomentosa above). Recognition of these two species is also justified from the perspective of the 
diagnostic species concept (Wheeler & Platnick 2000; Davis <Sr Nixon 1992), and the phenetic/taxononuc 
species Concept ( Judd ot al. 2007; Smith 1994) as evidenced by the data presented in Duncan and Brittain 
I '66). Finally, they must surely represent distinct evolutionary lineages (de Quieroz 2007; Wiley 1978)** 
evidenced by their reciprocal monophyly and difference in chromosome number (i.e., n - 12 in G. nana ,« 
= 24 in G. tomentosa). 


Given that Gaylussacia tomentosa and G. nana are largely sympatric, and very frequently co-occurnng, 
it is not surprising that we found a few morphologically intermediate specimens—likely hybrids—that ah 
exhibited polymorphic ITS sequences (possibly as a result of their possessing ITS sequences representatu* 
af the parental species). The frequency of such plants appears to be quite low (i.e., 7 specimens out o 



200 observed); most individuals (both m the field and herbarium) are easily determined. More research IS 
necessary on the breeding system of these plants, and the ploidy levels of both the parental species andtnt 
putative hybrids needs to be assessed. Despite the occurrence of these putative hybrids, we believe that o 
nana and G. tomentosa are largely reproductively isolated (possibly by the difference in chromosome numbe: 


thus fitting the biological species concept (Mayr 1969). 

The results presented here provide strong support for the treatments of Small (1933), Camp (1935,19™ 
Wood (1961), Duncan and Brittain (1966) and Luteyn et al. (1996), all of whom distinguished G. nana 
G. tomentosa at the species level. Recognition at the varietal level, as proposed or adopted by Gray (1»™ 
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Chapman (1889). Harper (WOo). Radford et 


L 4 


L964), Sleumer (1967), and Floyd (2002) is cc 



ippropriate. A mono 



Gavlussacia tomentosc 1 is present in the nrlTS and combined strict consensus 


rees albeit without BS support, although the morphological analysis does not resolve G. tomentosa into a 
lade in all trees (see Fig. 3). Twelve morphological characters distinguish G. tomentosa from G. nana (Table 




#2-9, 11-13). Five of the characters also separate G. tomentosa from G. JronJosa (chars. #2, 3, 4, 6 
and 7); these vary infraspecifkally to a greater extent within G. tomentosa and G. frondosa than in G. nana, 
except branch hair length (char. #3, Table 3). Finally, the difference in ploidy distinguishes G. tomentosa, 
which is the only tetraploid in the subsection (Love 1976; Luteyn et al. 1996), providing a distinctive apo- 
morphic feature for this taxon. 

The taxonomy of Gaylussacia frondosa is more problematic. Gaylussacia frondosa has no support as 
a monophyletic group in any of the strict consensus trees, yet the monophyly of this taxon is supported 
in some of the most parsimonious trees resulting from the combined analysis (Fig. 4). The two character 
states that could be considered diagnostic for G. frondosa are the absence of adaxial leaf hairs except for 
midvein and margin (char. #4; Tables 2, 3) and medium-length abaxial leaf hairs (#23). However, their 
pattern of variation is homoplasious. Thus, G. frondosa lacks clear morphological autapomorphies, which 
likely contributes to its lack of support in the nrlTS and morphological analyses. Although placed as a close 
relative of G. tomentosa in the combined analyses, we note that these two taxa can be clearly differentiated 
by chromosome number, n = 12 in G. frondosa and n = 24 in G. tomentosa (Love, 1976; Luteyn et al., 1996). 
In addition, G. frondosa has a more northern distribution (to NH and MA; Luteyn et al. 1996; Camp 1935) 
than does G. tomentosa, so in large parts of their ranges they are allopatric. It is commonly reported that 
t /wndosa can often reach heights of 2-3 m (Luteyn et al. 1996; Camp 1941), while G. tomentosa seldom 
teaches heights above 1.5 m (Luteyn et al. 1996). Camp (1941) also mentioned the highly clonal habit of G. 

''ttutiiosa and G* nana, as opposed to the more spread out, less rhizomatous habit of G. jrondosa. 

Most putative Gaylussacia tomentosa x nana hybrids were not included in the nrlTS analyses due to their 
-niorphic sequences, however, Gajdeczka 43 was included (Figs. 3, 4) as its ITS sequence is fairly clear. 


specimen (as discussed under results) is morphologically intermediate between G. tomentosa and G 


wwle.g., it has intermediate-length hairs on the abaxial leaf surface, and it may have been slightly glaucous) 
and its nrlTS sequence shows a slight indication of polymorphism at the variable nucleotide positions (he., 
nail secondary peaks are present at the polymorphic sites). The ITS sequence of this specimen was coded 
using the major peaks at these sites, and thus it was placed in the G. tOfHcntosdfftffndosu clade (instead ol 
ae G * nfl ^ clade) in the ITS tree (Fig. 3B). Had these sites been coded as polymorphic its position in the 
°8 ram would have been unresolved (witFin the G. /rondosa complex). Its more isolated placement in the 
net consensus of the combined analysis (Fig. 4), i.e., sister to the remaining members of the G. tomentosa 
G ‘frondosa clade, results from the interplay between the ITS sequence (with characters arbitrarily coded 
5 those of G. tomentosa) and the presence of certain morphological characters of G. nana. e.g., shottei hairs 
the abaxial leaf surface. Gajdeczka 43 is placed in the G. tomentosa clade (Fig. 3) in the morphology-based 
sis due to the more numerous apomorphies of this species (and indeed, this specimen is phenetically 
°si similar to G. tomentosa). Another putative hybrid, Gajdeczka 21 , was included in the morphology-based 
ylogenetic analysis, and this plant, like Gajdeczka 43, was placed in the G. tomentosa clade, likely as a 
t of its expressing more synapomorphies of G. tomentosa than G. nana (as coded in our matrix). Of 
S ^ 1 placements are influenced by the necessarily somewhat arbitrary delimitations of some of the 
VP ological character states (see methods). The nrlTS sequence of Gajdeczka 21 allows us to confidently 
dpo hesLe that this specimen represents a hybrid individual, and its placement with G. tomentosa 

J 1 mor Phological strict consensus tree (Fig. 3A) does not cause us to doubt that it is very likely of 


Such placements of hybrid individuals are to be expected (McDade 1990, 1992). 


}!ussacia baccata and G. WtsiHQ. were considered outgroups in this study, and as such, their intraspe- 
^ -Homy was not assessed (Floyd 2002). In addition, their species limits have not been controversial, 
flub as a cladospecies in the total evidence trees (Fig. 4). Strangely, G. baccata is paraphyletic 

Morphology-based trees (forming the basal branches of the cladogram). 
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The nrlTS analysis resulted in some unexpected results, e.g.. Gaylussacia f toccata and G. Ursirut formeda 
clade (Fig. 3B), although the resolution of the analysis was too poor to shed light on subsectional relationships, 
likely as a result of insufficient phylogenctically informative characters. Floyd’s (2002) nrlTS topology wj 
similarly unexpected (i.e., sect. Daamerium was shown as polyphyletic), and her nrlTS data was incompatible 
with her trnL-trnF data (as assessed by a partition homogeneity test). Possible lineage sorting or hybridiza¬ 
tion events affecting the nrlTS region might be contributing to the observed irregularities (Floyd 2002), and 
it should be noted that Floyd coded indels as additional characters, while we did not. Understanding the 
discrepancies between nrlTS sequences and other data sets used (or hypothesizing phylogeny in Guylussuud 
requires additional study. 

Our morphological ingroup topography also shows discrepancies with the nrlTS and combined topi 
raphies, i.e., Gay/ussucici frondosa groups with G. nana (rather than G. tomentosa ) in the morphological trees 
Three character states are synapomorphies of the G. frondosa + G. mina clade in the morphological tree: habit 
(#2), short nan-glandular hairs on the stems (#3), and moderately dense hairs on the abaxial leaf surface (07). 
Floyd’s (2002) morphological analysis places a G. tomentosa + G. nana clade (BS = 73%) sister to G. frondosa; 
G. tomentosa and G. nana also group together in Floyd's trnL-trnF topology. Her analysis, however, did not 
include any morphological character states that were shared by G. nana and G. frondosa, but also differed 
horn those of G. tomentosa. Our characters #2 and 3 are not treated in Floyd’s study, while character #7 is 
coded similarly for G. frondosa (but the state is shared with G. nana), and we note that it is not unexpected 
in analyses involving different evolutionary units to employ different morphological characters, or stair 
delimitations. Also in her analysis, several character states relating to hair density (glabrous twigs, adaxtal 
leal surfaces and leaf margins, few hairs on abaxial leaf surface) are unique to G. frondosa, and thus separai 
ii Irom G. fomcn/osu and G. nana. Twig hair density was measured in our study (and G. frondosa does have 
the least dense hairs), but this character was excluded from our analyses because it follows a similar pattern 
of variation as other hair density characters. It is possible that glabrous to sparsely pubescent twigs could 
represent a synapomorphy for the populations of G. frondosa. 

Finally, the morphological iBS = 85%) and combined (BS - 62%) data support the recognition of a 
monophyletic Gaylussacia frondosa complex (i.e., subsect. Frondosae) as it is commonly treated (Camp 1941 
Sletimer 1967; Floyd 2002). In the majority rule consensus tree of the nrlTS analysis (not shown), a G. W 
COta + G. ursina clade is sister to a G. tomentosa + G. frondosa clade in 52% of trees (with G. nana sister to all 
four species). In Floyd (2002) a G. tomentosa + G. frondosa + G. nana clade was evident in most analyses 

In conclusion, within Gaylassada sect. Decamerium subsect. Frondosae, it is clear that the largely sym- 
patric G. narnt and G. tomentosa each represent distinct species (regardless of the species concept a 
while the evidence for the distinctiveness of G. frondosa is more problematic—although our analyses pr 
some preliminary evidence that it, too, should be recognized as specifically distinct (Fig. 4). An identincatte 
key for these species is available in Luteyn et al. (1996). More work, however, needs to be done todetermir 
relationships within the complex, as the analyses are incongruent in the placement of G. frondosa , the rep¬ 
resentatives ol which are more closely related to G. nana in the morphological analysis and are more closet 
related to G. tomentosa in the nrlTS analysis. We also recommend additional investigations on the pk'idy 



level of members of subsect. FtXffldosete, especially the putative hybrids. 
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